New jamming scenario: From marginal jamming to deep jamming 
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We study properties of jammed packings of frictionless spheres over a wide range of volume 
fractions. There exists a crossover volume fraction which separates deeply jammed solids from 
marginally jammed solids. In deeply jammed solids, all the scalings presented in marginally jammed 
solids are replaced with remarkably different ones with potential independent exponents. Cor- 
respondingly, there are structural changes in the pair distribution function associated with the 
crossover. The normal modes of vibration of deeply jammed solids also exhibit some anomalies, 
e.g. strengthened quasi-localization and absence of Debye-like density of states at low frequencies. 
Deeply jammed systems may thus be cataloged to a new class of amorphous solids. 
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Colloidal suspensions and granular materials jam into 
amorphous solids when they are so compact that there 
is no room for constituent particles to move freely. As a 
simplified model to study the formation of this rigidity, 
a packing of frictionless soft spheres undergoes the jam- 
ming transition labeled "J" at a critical volume fraction 
4> c [lrl5|- Point J exhibits unusual criticality in the pres- 
ence of diverging length scales 043]) vanishing length 
scales in the pair distribution function g(r) [14l - ll6l |. and 
critical scalings of a variety of quantities which depend 
on inter-particle potential 0, [2 H, El, [13 ■ The nor- 
mal modes of vibration in marginally jammed solids near 
Point J also possess some anomalous properties. For in- 
stance, low-frequency modes are quasi-localized, anhar- 
monic, and poor in energy conduction [l3|, HHGjI. Some 
special vibrational features have been recently observed 
in experiments as well [20h23| . 

Although great efforts have been made to understand 
the anomalous properties of marginally jammed solids, 
it has not yet been questioned whether these properties 
would persist when the volume fraction keeps increasing 
away from Point J. In this letter, we extend the study 
of jammed solids to high volume fractions. The criti- 
cal scalings well-known in marginally jammed solids no 
longer hold in dense systems. Surprisingly, the scalings, 
structure, and vibrational properties of jammed solids 
undergo significant changes at approximately the same 
crossover volume fraction Moreover, properties of 
normal modes of vibration in deeply jammed solids at 
<t> > 4>d are beyond our understanding of both normal 
solids like crystals and marginally jammed solids. 

The three-dimensional systems studied here consist 
of a 50 : 50 binary mixture of N = 1000 frictionless 
spheres interacting via a spring-like repulsion: V (r^ ) = 
e (1 — rij/<Tij) a I a when the inter-particle separation of 
particles i and j, is less than the sum of their radii, 
<?ij = (crj + Oj)/2, and otherwise. Periodic boundary 
conditions are applied in all directions. The diameter ra- 
tio of the two species is 1.4 to avoid crystallization [24j . 
We set the small particle diameter a, particle mass m, 
and e to be one. The frequency is thus in the units of 
y/ e/ma 2 . 



It has been shown that 4> c is protocol dependent 
[25|,[26j]. However, critical scalings of <f> — <f> c (or equiva- 
lently pressure p) observed in marginally jammed solids 
are independent on the value of (j) c [2(|. We thus gen- 
erate jammed states at desired pressures to control their 
distance from Point J. We quench a random state at in- 
finitely high temperature and at an initial volume frac- 
tion close to that at the desired pressure to the local 
energy minimum by applying L-BFGS energy minimiza- 
tion method [27|- If the pressure of the jammed state 
is higher (lower) than the desired value, we decrease (in- 
crease) the volume fraction by a small amount 8<f> and ap- 
ply L-BFGS again. We apply this dilation-compression 
process successively until a jammed packing of spheres 
at desired pressure p is obtained. Note that 8<f> needs 
to decrease once the dilation and compression switch in 
order to be close enough to the desired pressure. At each 
pressure, we generate 1000 distinct jammed states from 
independent random states at infinite temperature and 
take averages over them. 

Figure [1] shows the potential energy per particle, V, 
bulk modulus, B, shear modulus, G, and coordination 
number, z versus the volume fraction <f) for both har- 
monic (a — 2) and Hertzian (a = 5/2) systems [28|]. In 
marginally jammed solids (<j) < <f>d), our data are fitted 
well with the following well-known scalings [5j: 
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where C w 0.6446 for a = 2 and 0.6443 for a = 5/2, 
and z c — 2d is the isostatic value at Point J with d the 
dimension of space. Most of these scalings depend on the 
inter-particle potential, which is one of the particularities 
of the jamming transition at Point J. Note that the fits to 
the bulk modulus using Eq. (Q]) slightly deviate from the 
data, as shown in Fig. [lib) . It is because that the bulk 
modulus in Eq. |T]) is derived from the widely accepted 
relation p ^ ((f) — <j> c ) a ~ 1 , but the actual exponent is 
slightly larger than a — 1 [29| . 

The scalings described in Eq. |T]) break down in the 
vicinity of a crossover volume fraction <pd independent 
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0.019 ± 0.004 
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TABLE I: Parameters S(<j>d) and us in Eq. @ for both harmonic and Hertzain systems. 
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FIG. 1: (color online) (a) Potential energy per particle V, 
(b) shear modulus G (empty symbols), bulk modulus B (solid 
symbols) , and (c) excess coordination number above isostatic- 
ity z — z c versus (f> — (j> c of jammed systems with harmonic 
(a — 2, circles) and Hertzian (a = 5/2, squares) inter- 
actions. The red dashed (harmonic) and blue dot-dashed 
(Hertzian) lines are fits for marginally jammed solids accord- 
ing to Eq. (0, while the solid curves are fits for deeply jammed 
solids according to Eq. ([2]). The insets to (a) and (c) show 
the half height width of the potential energy distribution at 
fixed pressure, SV, and the relative spatial variation of the 
coordination number, Sz/z versus (fi — fa, respectively, with 
the lines to guide the eye. 

on the potentials studied here. When <j> > fa, Eq. ([T]) 
turns into a unified scaling with cf> — (fid'. 

s-s{fa)~{fa-fays, (2) 

where fa ~ 1.180 for both potentials, and S denotes V, 
B, G, and z. The parameters S{4>d) and vs are shown 
in Table HI We can tell that the exponents vs do not 
depend on the inter-particle potential. 




FIG. 2: (color online) Pair distribution function of large parti- 
cles, g(r) for harmonic systems. The black solid, red dotted, 
green dot-dashed, and blue dashed lines are at <p = 0.688, 
1.107, 1.300, and 2.323. The inset shows g(r) at r = cr+, 
g(u^) (squares) and the first peak, gf^ (circles) versus 4>~4> c , 
with the lines to guide the eye. 



At fixed pressure, the distribution of the potential en- 
ergy of distinct jammed states is Gaussian. The half 
height width of the Gaussian distribution, SV is plot- 
ted versus fa- fa in the inset to Fig. [TJa). Interestingly, 
SV behaves different volume fraction dependence on both 
sides of fa: 8V ~ {fa- fa) a when < fa, and roughly a 
constant otherwise. 

In marginally jammed solids, spatial fluctuations are 
significant, e.g. the distribution of local coordination 
number is broad [13, IHj]. In the inset to Fig. QJc) we 
plot the relative spatial fluctuation of the coordination 
number, Sz/z versus <f> — fa. When 4> < (fid, Sz/z de- 
creases with increasing (fi and drops to the minimum at 
fa, while it increases with (fi when (fi > fa. Fluctuations 
in deeply jammed solids are thus not to be neglected. 

Marginally jammed solids undergo some structural 
changes in the pair distribution function g(r) during 
the jamming transition 

[T3-E1- The first peak of g{r), 
j™ x behaves a power law divergence approaching Point 
J. Meanwhile, the second peak splits into subpeaks at 
r = V3 and 2 in unit of particle diameter, which become 
discontinuous at Point J. Another anomalous feature is 
the discontinuity at r = 1 above Point J, implying that 
pairs of particles that are just in contact with each other 
are more than those that are almost in touch [14| . 

There are also structural signatures in g(r) associated 
with the crossover at fa, as shown in Fig.[2]of g(r) of large 
particles with harmonic interactions 30]. When (fi < <fid, 
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FIG. 3: (color online) Density of states, D(lj) of jammed 
solids with harmonic interactions. From the left to the right, 
the solid curves are measured at <j> = 0.6452 (red), 0.6508 
(orange), 0.688 (green), 1.107 (blue), 2.323 (magenta), and 
5.125 (maroon). The dashed line shows the Debye behavior 
with a slope of 2. The inset shows the exponent 7 of the 
power law fit, D(oj) ~ ui"' to the low frequency part of D(lj) 
versus <f> — </>c, with the line guide to eye. The error bars are 
smaller than the size of the symbols. 



there is only one peak in g(r) at r < <jl, where er^ is 
the large particle diameter. At <f>d 1 a second peak starts 
to emerge on the left-hand side of r = <jl , which locates 
at almost the same r in unit of ul when the volume 
fraction increases. Moreover, the inset to Fig. [2] shows 
that both the first peak, g™ 11 * and the right-hand side 
of r = ax, g(cr^) reach their minima at <pd- These three 
robust changes in g(r) structurally distinguish deeply and 
marginally jammed solids. 

A special feature of marginally jammed solids is the 
presence of a plateau in the density of vibrational states 
D(lu) j\. The plateau extends to cj — at Point J, im- 
plying the existence of excess soft modes at the jamming 
transition. It has been argued that D(uj) deviates from 
the Debye law {D(uj) ~ uj d ^ 1 ) especially close to Point J 
Q. However, for finite size systems, low- frequency modes 
are sparse and discrete in frequency, leading to difficulties 
in determining the exact behavior of the low-frequency 
D(uj). Furthermore, due to finite size effects, plane- 
waves hybridize with anomalous modes at low frequen- 
cies @, ITU . It is still a mystery if plane- wave- like modes 
could eventually dehybridize with anomalous modes and 
recover the Debye-like D(lj) at the low-frequency end in 
the thermodynamic limit. 

We diagonalize the Hessian matrix using ARPACK [3l[ 
to obtain the normal modes of vibration. Figure [3] shows 
the density of states D(co) for both marginally and deeply 
jammed solids with harmonic interactions [30j. D(ui) at 
the low frequency end can be approximately fitted with 
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The inset to Fig. [3] does not support the 
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FIG. 4: (color online) (a) Participation ratio, -P(i^) of jammed 
solids with harmonic interactions at <f> = 0.688 (black cir- 
cles), 1.107 (red squares), and 5.125 (blue triangles), (b) 
Normalized correlation of the polarization vectors, Cn{t) of 
the lowest frequency mode at <f> = 0.6508 (circles), 0.688 
(squares), 1.107 (left triangles), 2.323 (crosses), and 5.125 
(diamonds). The solid lines are fits according to Cjv(r) = 
Coexp(— r/£) + A. The inset to (b) shows £ versus <j> — <f> c , 
with the line to guide the eye. 



close to the jamming transition at Point J. It increases 
with increasing 0, but does not stop at 7 = 2. Inter- 
estingly, 7 reaches its maximum (~ 4) at the crossover 
volume fraction <pd- In deeply jammed solids, 7 decreases 
with increasing cj>. Debye behavior is completely absent 
in deeply jammed solids studied here. 

The anomalous density of states comes along with 
quasi-localization at low frequencies. Figure SJa) shows 
the participation ratio of normal modes [30(, defined as 



where is the polarization vector of 



Debye picture, i.e. 7 = 2. Our data indicate that 7 < 2 



particle i in the mode at oj, and the sums are over all 
particles. In both marginally and deeply jammed solids, 
there are always quasi-localized modes with low P(u>) at 
low frequencies. In marginally jammed solids, the quasi- 
localization is weaker with increasing 4>, expressed in the 
decrease of the number of quasi-localized modes and the 
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increase of the minimum P(u>). Figure |4j a) shows that 
near (fid there are quite a few low-frequency modes with 
large participation ratio forming a bump in P(uj), indi- 
cating the increase of the plane- wave components. When 
<fi > 4>d, however, quasi-localization is strengthened with 
increasing (fi. 

To quantify how the quasi-localization varies with the 
volume fraction, we measure the correlation of polar- 
ization vectors in the lowest frequency mode: C(r) = 
(e^ mm • e^ min ), where (...) denotes the average over con- 
figurations and all pairs of particles i and j with a sep- 
aration of r. Figure HJb) shows that the normalized 
correlation Cn(t) = C(r)/C(0) can be well fitted with 
C N (r) = C exp(-r/£) + A, where C > and A < 0. 
The parameter £ characterizes the size of the quasi- 
localized regions. The inset to Fig. 0|b) shows that £ is 
less than four particle sizes and reaches its maximum at 
the crossover volume fraction (f>d, so the lowest frequency 
mode at (fid is the least localized. 

All the results shown in this letter suggest that deeply 
jammed solids at <fi > (fid are new amorphous materials 
other than normal solids or marginally jammed solids. It 



is intriguing that most properties of jammed systems un- 
dergo remarkable changes at the same crossover volume 
fraction. The special features of deeply jammed solids 
must have significant impact on dynamics of thermal 
and sheared systems at very high volume fractions. We 
observe that the glass transition temperature increases 
with the volume fraction up to (fid and drops afterwards 
[32j , in consistent with a recent observation that increas- 
ing the density lowers the glass transition of ultrasoft 
colloids [IH, which is probably related to the correla- 
tion between the basin energy barrier height and quasi- 
localization [19j . We would expect similar volume frac- 
tion dependence of the yield stress. Based on this pic- 
ture, the jamming phase diagram [l|, [l|| may need cor- 
rections at high volume fractions for systems with soft 
potentials. Deeply jammed solids are accessible in ex- 
periments of core-softened colloids [HI . We hope that 
the present work could open a route to explore this new 
kind of amorphous materials. 
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